Effects of Gauss-Bonnet term on the final fate of gravitational collapse 
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We obtain a general spherically symmetric solution of a null dust fluid in n(> 4)-dimensions in 
Gauss-Bonnet gravity. This solution is a generalization of the n-dimensional Vaidya-(anti)de Sitter 
solution in general relativity. For n = 4, the Gauss-Bonnet term in the action does not contribute to 
the field equations, so that the solution coincides with the Vaidya-(anti)de Sitter solution. Using the 
solution for n > 5 with a specific form of the mass function, we present a model for a gravitational 
collapse in which a null dust fluid radially injects into an initially flat and empty region. It is found 
that a naked singularity is inevitably formed and its properties are quite different between n — 5 
and n > 6. In the n > 6 case, a massless ingoing null naked singularity is formed, while in the n — 5 
case, a massive timelike naked singularity is formed, which does not appear in the general relativistic 
case. The strength of the naked singularities is weaker than that in the general relativistic case. 
These naked singularities can be globally naked when the null dust fluid is turned off after a finite 
time and the field settles into the empty asymptotically flat spacetime. 

PACS numbers: 04.20.Dw, 04.20.Jb, 04.50.+h 



I. INTRODUCTION 



The final fate of gravitational collapse is one of the 
most important unsolved problems in gravitation physics. 
In general relativity, singularity theorems were proved 
stating that spacetime singularities inevitably appear un- 
der general situations and physical energy conditions [l| . 
In this context, a cosmic censorship hypothesis (CCH) 
was proposed by Penrose, which asserts that singulari- 
ties formed in generic gravitational collapse of physical 
matter cannot be observed; in other words, there are no 
naked singularities formed in physical gravitational col- 
lapse 0, Q- There are two versions of this hypothesis. 
The weak CCH asserts that all singularities formed in 
gravitational collapse are hidden inside black holes and 
cannot be observed from infinity, which implies the fu- 
ture predictability of the spacetime outside the black-hole 
event horizons. On the other hand, the strong CCH as- 
serts that no singularities are observable. 

At the present time, however, the general proof of CCH 
is far from complete; on the contrary, many counterex- 
ample candidates have been proposed in general relativ- 
ity |j| . The generic occurrence of naked singularities has 
been shown in the spherical dust collapse represented by 
the Lemaitre-Tolman-Bondi (LTB) solution HH El 
and in the null dust collapse represented by the Vaidya 
solution jjj 0, fill Iu3 ] although the absence of pressure 
is not physically reasonable. For a perfect fluid with 
a reasonable equation of state p = kfi, where p and /i 
are the pressure and the energy density, respectively, 
the generic occurrence of naked singularities in spheri- 
cal collapse has been shown by numerical simulations for 
< k < 0.0105 O Q O- Of course, the assumption 
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of spherical symmetry is very strong. In a genuine coun- 
terexample of CCH, the solution must be stable against 
non-spherical perturbations. Stability analyses against 
non-spherical linear perturbations have been done for the 
marginally bound LTB solution and the self-similar 
Vaidya solution ^3] ■ It nas been shown numerically that 
the naked singularity formation in the LTB solution is 
marginally stable against odd-parity perturbations but 
is unstable against even-parity perturbations |l6l |. On 
the other hand, the naked singularity formation in the 
self-similar Vaidya solution has been shown to be stable 
against even-parity perturbations [T^ |. In the perfect- 
fluid case, stability of the naked singularity formation 
against non-spherical perturbations is still open. 

Singularity theorems imply that a stage with infinitely 
high curvature can be realized in the very final stage 
of gravitational collapse. The effects of quantum grav- 
ity will be unavoidable there. Until now, many quantum 
theories of gravity have been proposed. Among them, su- 
perstring theory is the most promising candidate, which 
predicts a higher-dimensional spacetime. If superstring 
theory is the correct quantum theory of gravity, the ef- 
fects of the extra dimensions will become significant near 
singularity formation. When the curvature radius around 
a center is comparable to the compactification radius of 
the extra dimensions, the collapse must be effectively 
higher-dimensional. 

However, the non-perturbative aspects of superstring 
theory are so far not understood completely although 
the progress in recent years has been remarkable. Given 
the present circumstances, taking string effects pertur- 
batively into classical gravity is one possible approach to 
study the quantum effects of gravity. The Gauss-Bonnet 
term in the Lagrangian is the higher curvature correc- 
tion to general relativity and naturally arises as the next 
leading order of the ev'-expansion of heterotic superstring 
theory, where a 1 is the inverse string tension |l8j . 



2 



Until now, the properties of static solutions in Gauss- 
Bonnet gravity have been intensively studied, particu- 
larl y fr om the viewpoint of black-hole thermodynam- 
ics |l9l I20L I2TL I22I l23 |. In comparison to this, gravita- 
tional collapse has been less investigated, and accordingly 
the effects of the Gauss-Bonnet term on the dynamical 
formation process of such a higher-dimensional black hole 
have not been clarified yet. In this paper, considering 
the spherically symmetric gravitational collapse of a null 
dust fluid with the n(> 4)-dimensional action including 
the Gauss-Bonnet term for gravity, we investigate its ef- 
fects on the final fate of gravitational collapse. 

The organization of this paper is as follows. In Section 
II, we introduce our model and obtain an exact solu- 
tion which represents the gravitational collapse of a null 
dust fluid. In Section III, considering the situation in 
which a null dust fluid radially injects into an initially 
flat and empty region, we show that a naked singularity 
is inevitably formed. In Section IV, the strength of the 
naked singularity is investigated. Section V is devoted 
to discussion and conclusions. Throughout this paper we 
use units such that c = 1 . As for notation we follow [24| . 
The Greek indices run fj, = 0, 1, • • • ,n — 1. 



II. MODEL AND SOLUTION 

We begin with the following n-dimensional (n > 4) 
action: 



where 



■r-ir(R — 2A + aLas) 



r, (2.1) 



where R and A are the n-dimcnsional Ricci scalar and 
the cosmological constant, respectively. K n is defined 
by K n = ^/87rG rl , where G n is the n-dimensional grav- 
itational constant. The Gauss-Bonnet term Lqb is the 
combination of the Ricci scalar, Ricci tensor R^, and 
Ricmann tensor i? M v as 



Lgb = R 2 — 4i?, t „i? M1/ 



73 Wvpa 
±ij ixvpcr ±tj 



(2.2) 



In 4-dimcnsional spacctime, the Gauss-Bonnet term docs 
not contribute to the field equations, a is the coupling 
constant of the Gauss-Bonnet term. This type of action 
is derived in the low-energy limit of heterotic superstring 
theory [l8f. In that case, a is regarded as the inverse 
string tension and positive definite. Therefore, only the 
case where a > is considered in this paper. We con- 
sider a null dust fluid as a matter field, whose action is 
represented by Scatter in Eq. (|2.1|l . 

The gravitational equation of the action (|2.1|l is 



AS", 



kIT", 



(2.3) 



Hfii> — 2 



R, al3l R 



(2.4) 



(2.5) 



= holds for n = 4. The energy- momentum tensor 
of a null dust fluid is 



(2.6) 



where p is the non-zero energy density and is a null 
vector. 

We consider the general spherically symmetric spacc- 
time with the line element 

ds 2 = -f(v, r)e- 25iv ' r Uv 2 + 2 e - s dvdr + S(v, rfdn 2 n _ 2l 

(2.7) 

where dil 2 l _ 2 is the line element of the (n— 2)-dimensional 
unit sphere, v is the advanced time coordinate and hence 
a curve v = const, denotes a radial ingoing null geodesic. 



The normalization of 1^ is such that 1^ 



If S is constant, the (r, v) component of the field equa- 
tion H2.3fl gives the contradiction p = 0. Thus, we here- 
after adopt the gauge where S = r without loss of gener- 
ality. 

The (v,r) component of the field equation l|2.3(l gives 
5' = or / = I + r 2 /[2a(n — 3)(n — 4)], where the prime 
denotes the derivative with respect to r. In the latter 
case, which is valid only for n > 5, the (r, v) component 
of Eq. (|2.3() gives a contradiction p = 0, so that 5 = 5(v) 
is concluded. We can set i5 = without loss of generality 
by rescaling the time coordinate. 

Then, / is obtained by solving only the (v,v) or (r,r) 
component of Eq. I|2.3(l . which is written as 

rf'-(n-3)(l-f) + (n-l)Ar 2 

+ £(1 - /) [2rf - (n - 5)(1 - /)] = 0, (2.8) 

where a = (n - 3)(n - i)a and A = 2/[(n - l)(n - 2)]A. 
This equation is integrated to give the general solution 

as 



/ = 1 + 




(2.9) 



for n > 5, where m(v) is an arbitrary function of v. For 
n = 4, the solution is 



/ = 1- 



m(v) 



Ar 2 , 



(2.10) 



which is the Vaidya-(anti)de Sitter solution. In order 
for the term in the square root in Eq. (|2.9() to be non- 
negative, 



m > — I A + — 1 

4a , 



(2.11) 
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must be satisfied. There are two families of solutions 
which correspond to the sign in front of the square root in 
Eq. 1)2.9(1 . We call the family which has the minus (plus) 
sign the minus-branch (plus-branch) solution. From the 
(r, v) component of Eq. I|2.3|l . we obtain the energy den- 
sity of the null dust fluid as 



P 



2k 2 



(2.12) 



for both branches, where the dot denotes the derivative 
with respect to v. In order for the energy density to be 
non-negative, m > must be satisfied. 

In the general relativistic limit a — > 0, the minus- 
branch solution is reduced to 



/ = 1 



m(v) 



Ar 2 



(2.13) 



which is the ra-dimcnsional Vaidya-(anti)dc Sitter solu- 
tion, while there is no such limit for the plus-branch 
solution. In the static case m = 0, the solution ((2.9(1 
is reduced to the solution independently discovered by 
Boulware and Deser and Wheeler |2|j ■ 



III. NAKED SINGULARITY FORMATION 

In this and the next sections, we study the gravita- 
tional collapse of a null dust fluid in Gauss-Bonnet grav- 
ity and compare it with that in general relativity by use 
of the solution obtained in the previous section. We con- 
sider the minus-branch solution for n > 5 in order to 
compare with the general relativistic case. For simplicity, 
we do not consider a cosmological constant, i.e., A = 0. 
In this case, Eq. ((2.9|) is reduced to 




(3.1) 



Hereafter we call the solution the Gauss-Bonnet- Vaidya 
(GB-Vaidya) solution. The special case in which m is a 
constant we call the GB-Schwarzschild solution, of which 
the global structure is presented in j2^|. We consider 
the situation in which a null dust fluid radially injects at 

V = into an initially Minkowski region (m(v) = for 

V < 0) (Fig.^). The form of m(v) is assumed to be 



rriQV 



n — 3 



(3.2) 



where m is a positive constant. In this case, the so- 
lution ((3.1|> is reduced to the n-dimensional self-similar 
Vaidya solution in the general relativistic limit a — > 

o mm. 

From Eq. ((2.12(1 . it is seen that there is a central sin- 
gularity at r = for v > 0. The point of v = r = 
is more subtle but will be shown to be singular as well 
below. We will study the nature of the singularity. 

We have 



dr 
dv 



f_ 

2 



(3.3) 



along a future-directed outgoing radial null geodesic, so 
that the region with / < is the trapped region. A curve 
/ = represents the trapping horizon [23], i.e., the orbit 
of the apparent horizon, which is 



n — 3 



m(v) = ni{)V 
Along the trapping horizon, 



ds 2 



2(n — 3)mow™ 



a(n - 5)r"- 6 + (n - 3)r" 



-dv 2 



(3.4) 



(3.5) 



is satisfied and hence it is spacelike for v > and r > 0, 
i.e., it is a future outer trapping horizon, which is a local 
definition of black hole [23. From Eq. 1(3.4(1. only the 
point v = r = may be a naked singularity in the case 
of n > 6, while the central singularity may be naked for 
< v < i>ah in the case of n = 5, where uah is defined 
by 



m v AH 



(3.6) 



In order to determine whether or not the singular- 
ity is naked, we investigate the future-directed outgo- 
ing geodesies emanating from the singularity. It is shown 
that if a future-directed radial null geodesic does not em- 
anate from the singularity, then a future-directed causal 
(excluding radial null) geodesic does not also. (See Ap- 
pendix A for the proof.) So we consider here only the 
future-directed outgoing radial null geodesies. 

In order to show the existence of a null geodesic em- 
anating from the singularity, we adopt the fixed-point 
method 0, 0, We introduce a new coordinate 

■d = r/(v — vo) and then the null geodesic equation 13. 3|) 
becomes 



1 



dt) 



dv (v — vo) 



(d - i]) = i]^(v,d), 



(3.7) 



where we have introduced a parameter < 77 < 00 and 
i$ is given by 



(v - v ) 2 t} 2 




The constant vo is zero for n > 6, while it satisfies < 
^o 5: ^ah for n — 5. If we choose the parameter r\ to 
be 77 = rjo = (1 — \fmQVoj5i)/2 for n > 5, "J is at least 
C 1 in v > VQ^d > 0. Then we can apply the contraction 
mapping principle to Eq. ((3.7(1 to find that there exists 
the solution satisfying &(vq) = rj , and moreover that it 
is the unique solution of Eq. ((3.7(1 which is continuous at 
v = v Q . (See 0,0,111 for the proof.) 

Now we have shown that there exists a future-directed 
outgoing radial null geodesic which behaves as 



2r 



(3.9) 
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near v = r = for n > 6, while as 

v - v o + y 2 . 2/ . r (3.10) 

near r = and v = vq with < vq < «ah for n = 5. 
Along the null geodesic emanating from the singularity 
at u = r = for n > 5, the Kretschmann invariant 
A' = i? M „p CT i? A " /po , which is represented by / as 

(3.11) 

diverges for r — > as 

A = 0(l/r 2 ). (3.12) 

Along the null geodesic emanating from the singularity 
at r = and < v < «ah for tt, = 5, A diverges as 

A = 0(l/r 4 ), (3.13) 

and hence they are actually singular null geodesies. We 
have now shown that at least a locally naked singularity 
is formed, and consequently the strong version of CCH 
is violated. 

Next we consider the structure of the naked singularity. 
First we consider the case of n > 6. We expand the 
trapping horizon lj3.4J) around r = as 

v = mo 1/(n ~ 3 V"~ 5 +r"- 3 ) 1/( "~ 3) , (3-14) 
~ (5/m ) 1 /("- 3 )r 1 - 2 /(™- 3 ). (3.15) 

Thus, there exists a spacetime region hi which is both 
the past of the trapping horizon and the future of the 
future-directed outgoing radial null geodesic 7 which be- 
haves as Eq. (|3.9|l near v = r = 0. Because the trapping 
horizon is spacclike for v > and r > 0, the past-directed 
ingoing radial null geodesic C emanating from an event 
in hi never crosses the trapping horizon. Also £ never 
crosses 7 except for v = r = because they are both 
future-directed outgoing radial null geodesies. Conse- 
quently, f inevitably reaches the singularity at v = r = 0. 
Since U is an open set, it is concluded that there ex- 
ists an infinite number of future-directed outgoing ra- 
dial null geodesies emanating from the singularity at 
v = r = 0. Such geodesies should correspond to the 
solution of Eq. (|3.7() with 77 = or 77 = 00. On the other 
hand, the future-directed ingoing radial null geodesic ter- 
minating at v = r = is only v = 0, thus we conclude 
that the naked singularity at v = r = has the ingoing- 
null structure. 

In the case of n = 5, the singularity at r = and 
v = vq with < vq < «ah is both an endpoint of the 
future-directed ingoing radial null geodesic and the initial 
point of the future-directed outgoing radial null geodesic, 
therefore we conclude that the naked singularity has the 
timclike structure. 




r 



v=0 

FIG. 1: Gravitational collapse of a null dust fluid in (A) 5- 
dimension with an event horizon (vf > hah), (B) 5-dimension 
without an event horizon (vf < «ah), and (C) n(> 6)- 
dimension. A singularity is formed at v = r — and develops. 
Here the Minkowski spacetime for v < is joined to the GB- 
Schwarzschild spacetime for v > V{ by way of the GB-Vaidya 
spacetime. 



On the other hand, it cannot be concluded from the 
present analysis that there exist no null geodesies ema- 
nating from the singularity at r = and v = i>ah be- 
cause null geodesies which correspond to the solution of 
Eq. (|3.7|) with 77 = or 77 = 00 might exist. More de- 
tailed analyses are needed in order to clarify the nature 
of that singularity at r — and v — i>ah • If there is a null 
geodesic emanating from the point r = and v = uah, 
it is concluded by the same discussion as that in the case 
of n > 6 that the singularity has an ingoing null portion. 

We can consider the situation in which the null dust 
fluid is turned off at a finite time v = Vf > 0, where- 
upon the field settles into the GB-Schwarzschild space- 
time with m = niQv" ~ 3 , which is asymptotically flat. 
The Minkowski spacetime for v < is joined to the GB- 
Schwarzschild spacetime for v > V{ by way of the GB- 
Vaidya solution (|3.1|) (Fig. . Then the singularity is al- 
ways globally naked when n = 5 with V{ < wah since no 
horizon is formed. When n > 6 or 11 = 5 with Vf > v ah, 
the singularity can also be globally naked. If we set Vt 
to a sufficiently small value and a sufficiently close value 
to v ah for n > 6 and n = 5, respectively, a null ray ema- 
nating from the singularity reaches the surface v = V{ in 



5 



the untrappcd region, and consequently it can escape to 
infinity. In that case, the singularity is globally naked, 
i.e., the weak version of CCH is violated. In order to 
determine the values of the parameters with which the 
singularity is globally naked, numerical integrations of 
the null geodesic equation 1)3.3(1 are required, in which 
the asymptotic solutions (|3.9[) and (|3.1U|) are used as the 
initial conditions. 

The structures of the singularity in the outer GB- 
Schwarzschild spacetime are as follows j^. The singu- 
larity is spacelike for n > 6 and for n = 5 with Vf > «ah- 
It is timelike for n = 5 with Vf < wah, while it is null for 
n = 5 with V{ = uah- Together with these results, the 
possible global structures of the gravitational collapse are 
drawn in Fig. [2] 

Next we consider the general relativistic case, i.e., the 
?i-dimensional Vaidya solution. In this case, the trapping 
horizon is represented by v = v^r, where 



(A) r=0 



1 



l/(n-3) ' 



(3.16) 



and hence only the point v = r = may be naked. We 
find the exact power-law solution of the null geodesic 
equation 13.3fl : 



vir, 



when 



g(vi) = m < 2 - «i + 2 = 



(3.17) 



(3.18) 



has at least one positive real root. This condition is re- 
alized when 



g(v ex ) < o 

holds, where w ox = V[( n ~ 2)too] 1 ^™ -3 ^ satisfies 
dg(v e , 



dvi 



0. 



The condition H3.19|l is written as 

(n - 3)"~ 3 



mo < 



)n-3 



(n-2)» 



(3.19) 



(3.20) 



(3.21) 



It is emphasized that Eq. I|3.17|l is not an asymptotic so- 
lution but an exact solution of the null geodesic equation 
(|3.3|1 . Along this null geodesic, the Kretschmann invari- 
ant diverges for r — > as 



K = 0(l/r 4 ), 



(3.22) 



and hence it is concluded that an ingoing-null naked sin- 
gularity is formed at v = r = 0. This result is consistent 
with the results in pi . 

Here we show that if Eq. I|3.21|l is satisfied then the 
null geodesic l|3.17|l can reach infinity and therefore the 




FIG. 2: Possible global structures of the globally naked sin- 
gularity formation. Zigzag lines represent the central singu- 
larities, corresponds to the future (past) null infinity. 
BEH and CH stand for the black-hole event horizon and the 
Cauchy horizon, respectively. (See Ref. for rigorous defi- 
nitions of the event horizon and the Cauchy horizon.) Naked 
singularity formation in general relativity is represented by 
(A). In 5-dimension with the Gauss-Bonnet term, the global 
structure is represented by either (C) or (D) ((E) or (F)) for 
Vf > (=)uah, while it is represented by (B) for V{ < vah- 
In n(> 6)-dimension with the Gauss-Bonnet term, the global 
structure is represented by (A). There are counterparts of (A), 
(C), and (D) which represent the locally naked singularity for- 
mation. 
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singularity is globally naked independent of the value of 
Vf. Because of 



^- = {n- 2)m <" 3 
avi 



1 



(3.23) 



and g(0) = 2 > 0, the null geodesic equation (|3.3|l rep- 
resents either one or two solutions when Eq. (|3.21l) is 
satisfied. We denote them as v — v a r and v = v^r 
with < v a < v cx < Ufc, where we have equality when 
g(v e x) = 0. From Eq. 1)3.16)1 . Vh > v cx > v a satisfies, so 
that the null geodesic v = v a r is always in the untrapped 
region and reaches the future null infinity independent 
of the value of V[. Thus, when Eq. (|3.21() is satisfied, 
the solution represents the formation of a globally naked 
singularity, i.e., the weak CCH is violated. 

In the case that mo does not satisfy the condition 
(|3.21() . more detailed analyses are needed to determine 
whether the singularity is naked or censored because it 
might be possible that there exist null geodesies which 
do not obey the power law near v = r = 0. 



IV. CURVATURE STRENGTH OF A NAKED 
SINGULARITY 

In the previous section, it was shown that a naked 
singularity is inevitably formed in Gauss-Bonnet gravity. 
In this section, we investigate the curvature strength of 
the naked singularity. We define 



(4.1) 



where = dx^/dX is the tangent vector of the future- 
directed outgoing radial null geodesic which emanates 
from the singularity and is parameterized by an affinc 
parameter A. We evaluate the strength of the naked sin- 
gularity by the dependence of ib on A near the singularity. 

In the 4-dimensional case, the strong curvature con- 
dition (SCC) |3(il ] and the limiting focusing condition 
(LFC) l31l are defined. Following the work of Clarke and 
Krolak |33|, consider a geodesic (7) with tangent vector 
and terminating at or emanating from a singularity, 
where A = 0. If limA^o A 2 ?/' > and lim^o Xip > 0, 
then SCC and LFC are satisfied along 7, respectively. 
However, these results have not been extended to the 
higher-dimensional case so far. 

From Eq. (|3.3|l . a relation 



2 



(4.2) 



is obtained. Using this, we obtain from Eq. 1)4. l|l that 



4, 



rf 2 



(4.3) 



We write the geodesic equation for the radial null 
geodesic: 



dX 



r\2 



(4.4) 



First, we consider the case in Gauss-Bonnet gravity. 
Along the null geodesies 1)3. 9fl and 13.1Uf) . we obtain 



lim/ = 1, 



= _n_3 ./2-mo 

r^0 J 2 



for n > 6 and 



lim / = 1 - 

r^O V a 

lim/ = —\—r- 

r^0 J V a 



(4.5) 
(4.6) 

(4.7) 
(4.8) 



for n = 5, respectively. Using Eqs. H4.5() - (|4.8() . we find 
that Eq. H4.4f> is reduced to the following form near the 
singularity: 



where q is defined by 



q={n- 3) 



2"- 5 m 



for n > 6, while by 



2y/m Q /c 



(1 - y/m vl/a) 2 
for n = 5. The solution of Eq. I|4.9() is 



r ~ -hx{CX + D), 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



where C and D are integration constants. Here, we set 
C = D = 1 without loss of generality by redefinition of 
the affinc parameter so that A = corresponds to the 
singularity r — 0. 

After straightforward calculations using Eqs. (|4.5() ~ 
(|4.8|) and (|4.12|1 . we obtain from Eq. (£OJ) that 



n - 2 



for n > 5. Thus, 



lim ib — 

a^o Y ln(A + 1) 



lim Xib = n — 2 



(4.13) 



(4.14) 



is satisfied for n > 5. 

Next, we consider the case in general relativity. In 
this case, we obtain from Eq. 1)4.4)1 that the null geodesic 
()3.17)) obeys the null geodesic equation 



where 



Afcr = _£ (fc r )a 

dX r 



(n — 3)rriQVi 



(4.15) 



(4.16) 
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The solution of Eq. (|4.15|) is 

r = (E\ + Fy/ (1 + S \ 



(4.17) 



where E and F are integration constants. We set E = 1 
and F = by redefinition of the affine parameter. After 
calculations in a manner similar to those in Gauss-Bonnet 
gravity, we obtain 



(n - 2)s 



(1 



) 2 A 2 ' 



(4.18) 



This result is consistent with the results in |2a.l2a|. 

Consequently, it is concluded by comparing Eq. (|4.13|) 
with Eq. H4.18fl that the strength of the naked singularity 
in Gauss-Bonnet gravity is weaker than that in general 
relativity. 



V. DISCUSSION AND CONCLUSIONS 

We have obtained an exact solution in Gauss-Bonnet 
gravity, which represents the spherically symmetric gravi- 
tational collapse of a null dust fluid in n(> 4)-dimensions. 
For n > 5, the solution is reduced to the n-dimensional 
Vaidya-(anti)de Sitter solution in the general relativistic 
limit. For n = 4, the Gauss-Bonnet term does not con- 
tribute to the field equations, so that the solution coin- 
cides with the Vaidya-(anti)dc Sitter solution. Applying 
the solution to the situation in which a null dust fluid ra- 
dially injects into an initially Minkowski region, we have 
investigated the effects of the Gauss-Bonnet term on the 
final fate of the gravitational collapse. 

We have assumed that the mass function has the form 
of m(v) = mov n ~ 3 . Then, it has been found that there al- 
ways exists a future-directed outgoing radial null geodesic 
emanating from the singularity in Gauss-Bonnet gravity, 
i.e., a naked singularity is inevitably formed. On the 
other hand, in the general relativistic case, there exists 
such a null geodesic only when mo takes a sufficiently 
small value. This result implies that the effects of the 
Gauss-Bonnet term on gravity worsen the situation from 
the viewpoint of CCH rather than prevent naked singu- 
larity formation. 

Furthermore, the Gauss-Bonnet term drastically 
changes the nature of the singularity and the whole pic- 
ture of gravitational collapse. The picture of the gravi- 
tational collapse for n = 5 is quite different from that for 
n > 6. For n > 6, as well as the general relativistic case 
for n > 4, a massless ingoing null naked singularity is 
formed. On the other hand, for the special case n = 5, a 
massive timelike naked singularity is formed. It is found 
from Eq. (|3.4|) that the formation of a massive timelike 
singularity in 5-dimcnsions is generic for the general mass 
function m(v) which satisfies m(0) = and in > 0. 

Here, we should mention the substantial work by Lake 
in 4-dimensional spherically symmetric spacetimes |33j . 
It has been shown under very generic situations that the 
massive singularities formed from regular initial data are 



censored [33j. In other words, if a naked singularity is 
formed from regular initial data, it is massless or with 
negative mass. This result was obtained without using 
the Einstein equations and may be extended into higher 
dimensions. One might think, therefore, that this result 
is inconsistent with ours in Gauss-Bonnet gravity with 
n = 5. However, the discussion and results in [33| | can- 
not apply directly to our system. In Lake's analysis, the 
Misncr-Sharp mass is adopted to evaluate the mass of 
the singularity. On the other hand, we adopted m(v), 
which is different from the higher-dimensional Misner- 
Sharp mass. As shown in Appendix B, m(v) is prefer- 
able in Gauss-Bonnet gravity to evaluate the mass of the 
singularity. 

Although naked singularities are inevitably formed in 
Gauss-Bonnet gravity, the Gauss-Bonnet term makes the 
strength of the naked singularity weaker than that in the 
general relativistic case. In association with this, there 
does exist a possible formulation of CCH, which asserts 
that the formation of weak naked singularities need not 
be ruled out [34||. In this sense, the Gauss-Bonnet term 
works well in the spirit of CCH. 

Lastly, several notes should be made. In this paper, 



we have assumed m(v) = mov T ' 



for simplicity. Anal- 



yses with a more general form of m(v) are needed in 
order to determine whether the nature of the singular- 
ities and the picture of gravitational collapse obtained 
in this paper are generic or not. Also, the higher-order 
correction terms to gravity might change the nature of 
singularity and the final fate of gravitational collapse. In 
addition, the assumption of spherical symmetry is very 
restricted. Very recently, the naked singularity formation 
in the self-similar Vaidya solution was found to be sta- 
ble against non-spherical linear perturbations with even 
parity |17||. Analyses in Gauss-Bonnet gravity beyond 
spherical symmetry must be addressed. These ideas are 
under investigation. 
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Note added: 

After completing this work we were informed that 
Kobayashi also obtained the solution (|2.9|) and its gener- 
alization to the spacetimc M w M 2 x K, n ~ 2 where K, n ~' 
is the (n — 2)-dimensional Einstein space |3 
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APPENDIX A: GEODESICS FROM 
SINGULARITIES 

In this appendix, we show that if a future- directed 
causal (excluding radial null) geodesic emanates from 
the central singularity, then a future- directed radial null 
geodesic emanates from the central singularity. The con- 
traposition of this theorem implies that it is sufficient 
to consider only the future-directed outgoing radial null 
geodesies in order to determine whether or not the sin- 
gularity is naked. The proof is similar to that in the 
four-dimensional case in p3q . 

In the GB-Vaidya spacctime 



ds 2 = -fdv 2 + 2dvdr 



l dCl 



ra-2) 



(Al) 



where / is given by Eq. I|3.1[) . the tangent to a causal 
geodesic satisfies 



emanating from the singularity in the trapped region, 
i.e., the singularity with / < is censored. 

Next we focus on the singularity with < v < vau, 
where / > holds. Now suppose that v = vcg(?") ex- 
tends back to a central singularity located at (r, v) = 
(0,v s ), where v s satisfies < v s < wah- There exists a 
portion of v = ucgO") which is in the untrapped region. 
Let p be any point on such a potion of v = vqg ( r ) an d to 
the future of the singularity. Applying inequality l|A6|) at 
p, we see that v = vrng (?*) through p crosses v = i>cg(t) 
from above and hence the points onv = vrng ( r ) prior to 
p must lie to the future of points on v = vcgM prior to 
p, in the sense f rng ij" ) > v CG{ r ) for r £ (0, r*), where r* 
corresponds to p. Thus, the radial null geodesies must ex- 
tend back to r = with v = vq satisfying v s < Vq < «ah, 
and so must emerge from the singularity. 

It is noted that i>ah = v s = vq = holds in the case of 
n > 6. 



/ 



/ dv \ 2 dv dr 1? 
— +2—— + — 



\d\ 



dX dX 



(A2) 



APPENDIX B: MASS OF NAKED SINGULARITY 



where A is an affine parameter, L is the conserved angular 
momentum, and e = 0, — 1 for null and timclikc geodesies, 
respectively. Then, at any point on such a geodesic, 



> 2 



dv dr 
dXdX 



(A3) 



with equality holding only for radial null geodesies. For 
the future-directed outgoing geodesies, this gives 



and 



dv 2 

* s 7 >0 



dv 2 „ 

^7 <0 



(A4) 



(A5) 



in the untrapped and trapped regions, respectively. 
Therefore, 



dv C G cfoRNG Q 



dr 



dr 



and 



dv C G dvn^G Q 



dr 



dr 



(A6) 



(A7) 



are satisfied in the untrapped and trapped regions, re- 
spectively, where the subscripts represent causal (ex- 
cluding radial null) geodesies and outgoing radial null 
geodesies, respectively. 

Let us consider the (r, i>)-plane. The singularity is lo- 
cated at r = for v > 0. First we consider the singular- 
ity with / < 0. By Eq. (|A7|) . the past-directed ingoing 
geodesies emanating from an event with r > in the 
trapped region cannot reach the singularity at r = 0. 
Therefore, there is no future-directed outgoing geodesic 



In this appendix, we show that the higher-dimensional 
Misncr-Sharp mass is not an appropriate quasi-local mass 
in Gauss-Bonnet gravity. 

Firstly, we give a definition of the higher-dimensional 
Misner-Sharp mass. We consider the n-dimensional 
spherically symmetric spacetime M f=a M 2 x S n ~ 2 with 
the general metric 



M „ = diag(# AB^r 2 ^), 



(Bl) 



where gAB is an arbitrary Lorentz metric on M. , r is a 
scalar function on M 2 with r — defining the boundary 
of A4 2 , and jij is the unit curvature metric on S n ~~ 2 . The 
higher-dimensional Misner-Sharp mass ttims is a scalar 
on M 2 defined by 



(n - 2)K- 2 
2kI ' 



\l-r. A r A ). 



(B2) 



There are several properties which a well-defined quasi- 
local mass should satisfy [3Hl38l| . One of them states that 
there is an invariant mass function in spherically symmet- 
ric spacetime, that any definition of the quasi-local mass 
should be reduced to in the special case with spherical 
symmetry. In particular, in the Schwarzschild spacetime 
of mass Msch, the invariant mass function should coin- 
cide with Mg c h because it is the unique vacuum solution 
in general relativity. In Gauss-Bonnet gravity, on the 
other hand, the unique spherically symmetric vacuum 
solution is the GB-Schwarzschild solution. Therefore, 
a well-defined quasi-local mass in Gauss-Bonnet gravity 
should be reduced to the mass for this solution in the 
spherically symmetric vacuum case. 

The metric of the GB-Schwarzschild solution is given 

by 



ds 1 



^fdv 2 + 2dvdr 



r 2 dn 2 n _ 2 



(B3) 
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with 

/=i+ S(W i+4fi ^)' (B4) 

where m is a constant. The rt-dimcnsional Schwarzschild 
solution is obtained in the limit of a — > as 

/ = l-^5. ( B5 ) 

of which mass is given by Ms c h = (n — 2) V n _2ra/(2K^), 
where V^_2 is the volume of an (n — 2)-dimcnsional 
unit sphere. On the other hand, the mass for the GB- 
Schwarzschild solution also coincides with Ms c h plj - 
Thus, a well-defined quasi-local mass should be reduced 
to Msch for the spherically symmetric vacuum case in 
Gauss-Bonnet gravity as well as in general relativity. 



From Eq. (|B2|) . we can show toms = M^ch for the n- 
dimensional Schwarzschild solution, while we obtain 

for the GB-Schwarzschild solution. Indeed, tums does not 
coincide with Ms c h for n > 5 for the latter case. Conse- 
quently, ttims is not an appropriate quasi-local mass in 
Gauss-Bonnet gravity. 

In this paper, 777(7;) in Eq. is adopted in order to 
evaluate the mass of the singularity. 777(71) is preferable 
rather than 777ms m Gauss-Bonnet gravity because it is 
reduced to 2nl a Mschl[{n — 2)V n -2\ for the static limit, 
i.e., in the GB-Schwarzschild solution. 
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